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ABSTRACT 

We discuss the phenomenological extraction of the leading j -plane singularity from 
singlet structure functions F s measured at small x. Using a saddle-point method we show 
that dlni^ s /dln^ is a suitable observable for this purpose in the region x < 10 -2 . As an 
application, we confront and distinguish in a model-independent way structure function 
parametrizations coming from two different QCD evolution equations, namely the Lipatov 
(or BFKL) equation and the Gribov-Lipatov-Altarelli-Parisi (or GLAP). Recent results at 
HERA are discussed in this framework. 



Recent results from HERA experiments (H^ and Zeus^) yield information on the 
singlet proton structure function F s for a new range in the Bjorken variables, namely 
10 -4 < x ^ 10 -2 , 8 < Q 2 < 60 GeV 2 . These experimental results, although of low 
statistics, have already provided promising theoretical information^ . 

The rise of F s as x becomes very small, has reinforced predictions of parametrizations 
based on perturbative QCD evolution equations, of a large increase in space-like gluon 
radiation. The QCD calculation in the leading (o^ln 1 ;)™ approximation^ appears to 
be relevant for this phenomenon, even though non-perturbative effects^ or gluon-gluon 
interactions^ may modify its quantitative predictions. 
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In this framework, the structure function is dominated by a singularity in the j-plane 
of its Mellin transform, namely 

F s (x) ~ x~ a (in^J 

£ dx x^F s (x) ~ F a (j) ~ {j _ 1 a)b+1 (1) 

with a > 1 being the end-point singularity, and b characterizing the type of singularity 
(square root b = — |, pole 6 = 0, etc.). An exact calculation of the Lipatov (BFKL) 
singularity at a fixed external momentum scale yields b = — ^ but this coefficient may be 
substantially different for structure functions^ . 

The Lipatov singularity is given by a = 1 + 41n2^a s and may reach 1.5. However, 
even though this perturbative approach looks promising, some ambiguities have still to 
be resolved. For instance, non-perturbative contributions are present because of the large 
diffusion in the gluon transverse momentum scales contributing to the structure function^ . 
More intringuingly a prediction (GRV parametrization^) based on the GLAP evolution 
equation^ has also been successful at HERA. Its success is surprising because it is based on 
extrapolating from a very small value of Q 2 = /i 2 by summing only the leading [a s lnQ 2 ] n 
contributions, and thus on a different picture of the j-plane singularity. Is its success a 
mere coincidence or is it due to a subtlety within the theory? This is the question we wish 
to settle. 

Our aim in this letter is to answer this question by proposing a model-independent 
study of the j— plane effective singularity of the structure functions. Our results are as 
follows: 

i) When using a saddle-point approximation valid at small x (x < 10 -2 ), the Mellin 
transform of eq.(l) reduces to a local correspondence, x > j{x) between F s and F s (j) . 

ii) Within the same approximation, the function dlni^ s (x)/dln-| depends only on the 
rightmost j— plane singularity of F s . This result is in agreement with a recent empirical 
observation^ . 
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iii) Applying this result to the BFKL or GLAP j — plane singularity, we show that 
while different analytically, they lead to an almost identical dlni^ s (x)/dln-| in the small x 
range of the HERA experiments, at fixed Q 2 . 

iv) The only way to distinguish the two types of singularities may be the Q 2 depen- 
dence of dlni^ s (x)/dln-^. 

Let us prove assertion i). The saddle-point approximation relating a structure function 
at small x and its Mellin transform near its rightmost singularity is correct for the Lipatov 
singularity! 11 ^. \y e wm use this method in a general way. Let us consider the inverse Mellin 
transform of relation (1) 

Fs(x) = / e *"/*+fdj (2) 

2z7T J C 

where / = \nF s (j) and C is a contour in the complex j -plane with Re j = jo at the right 
of all singularities of F s (j). 

For ml/a: large enough, we get using the saddle-point approximation 

lnF a ( a ;)~jln- + /-^ln(27r/"), (3) 

x 2 

where the saddle point j(x) is defined by the condition 

= ln- + /'(j) (4) 

x 

and the primes mean the derivative w.r.t. j at j = j, i.e. /' = ^ . Equation (4) gives 

a local relation x > j( x )- Moreover, /' is given by the dominant j-plane singularity of 

F s (j). As an example the singularity given in eq.(l) yields j = a + (b + l)/ln^-. The value 
of [ln-j:] 1 is a measure of the distance to the singularity at j = a. 

Next we prove assertion ii) about dlni^ s /dln-|. Using the saddle-point formula (3), we 
get by differentiating, 



Note that f^r = [/"] • Finally, using (4), the derivative reads 



d InF 



T^j-hf"- 1 ]', (5) 



d ln± 2 

X 

(this simple expression stems from the fact that the structure function is a Legendre 
transform for the conjugate variables j and ln-j:, see eq.(2)). Now we note that eq.(5) is 
dominated by the rightmost j-plane singularity of F s . Indeed, non-singular (or less singular) 
contributions appears only in the saddle-point correction term, namely , and they 

are completely negligible in the small x domain. 

Indeed, a similar observation has been made recently about a background-subtracted 
structure function^ . Our derivation, using the saddle-point method, provides an expla- 
nation for this empirical observation. The background has a singularity near j ' = 1, which 
has to be subtracted out in order to extract the main BFKL contribution as done in ref. 
[9]. 

Now let us consider the well-known j— plane singularity generated by the GLAP evo- 
lution equations^ . Define our notation by 

~ ( Vp — V- v+ — Vp \ 

f s(j) = ex p u +z h ex p v -i (? + ?) (J) 

\U+-U- U+-U- J 

1N f <j)Z . ^ ^ ... 

H (exp v+i - exp V-Q g(j ) 



v + - v F 

Vp — V- 
V+ — V- 



exp v+f^ * (a regular function at j > 0) , (6) 



where £ is the well known evolution parameter of the GLAP equation and (q + q) (resp. g) 
are the input quark singlet (resp. gluon) structure functions at £ = 0. Note that here we 
want to analyze the singularity structure of the GLAP evolution equation. Thus, we will 
apply our formalism to the case of input functions q } q } g which are regular for j > 0. 
The GLAP kernels read 

4 

v F = --\ 4tP(j + 1) + 4j E - 3 



2 24 (j 2 + j + 1 

// G = -12 + 7^1 + 11--^/+ V 



3 ' j(j 2 -1)(j+2) 



j 2 +j + 2 lG 8j 2 +j + 2 



G i(j + l)(i + 2) 3 i 0"2 _ i) 



^ { f f + v G ± - z/ G ) 2 + 8N f <f>$<f>% 1 , (7) 



where ?/> is the logarithmic derivative of the T function and the Euler constant. 

At this stage, some comments are in order. The full leading singularity at j = 1 in the 
GLAP equations is taken into account since all singular contributions at j = 1 (the essential 
singularity contained in e l> +^ and its polynomial prefactor) are present. This is due to the 
fact that up, cf>Q are non- singular at j ' = 1 whereas v+ = ^j-+regular function at j > 0. 

On the other hand, the whole singularity structure of the GLAP equation depends 
only on the parameter £ 

In Q 2 /A 2 



where Q is the scale of the deep-inelastic reaction, /i, the initial evolution scale, A, the QCD 
one-loop scale and Nf the number of flavors. Note that a phenomenological application^ 
may require evolution beneath the perturbative domain i.e. (ji < Q < Qo), where Q\ ~ 
GeV 2 . In our study, we shall use £ to parametrize the GLAP singularity including the 
non-perturbative evolution range. 

Now let us demonstrate assertion iii), i.e. the equivalence between the GLAP and 
BFKL singularities in the small x range at fixed Q 2 . In Fig.l, the predictions of dlmF/dln-j: 
are given for three typical values of £ and compared with the 2-parameter fit using the 
BFKL expression (1). The a and b parameters of eq.(l) are displayed as a function of £ in 
Fig. 2. Surprisingly enough, although the singularities are different analytically, both give 
similar results in the region of interest for HERA, 10 -4 < x < 10 -2 . We have performed 
a variety of checks on this claim, varying the fitting range, and switching on and off the 
non-singular terms in the GLAP kernels themselves. This results in a small spread in the 
values of a and b, see Fig. 2. For x > 10 -2 the regular terms can no longer be neglected, 



while for x < 10 4 the difference between the two parametrizations becomes noticeable, 
though still small. 

Now let us come to our final assertion iv) about the Q 2 -evolution of dm-F/dm^-. The 
BFKL singularity is not expected to vary significantly with Q 2 , even if we take into account 
the running of the coupling constant^ 11 ! . To be more specific, the overall normalization of 
the singlet structure functions is indeed expected to vary as \J Q 2 but its shape is expected 
to remain unchanged^ 12 ]. This implies that dlmF/dln 1 ; is independent of Q 2 . In contrast, a 
distinctive feature of the GLAP singularity is its Q 2 evolution. For the quantity dlmF/dln 1 ; 
this evolution is determined by the £ dependence of a and b, see Fig. 2. Note that a Q 2 
dependence compatible with the first-order GLAP equations has also been obtained in 
phenomenological studies of parton distributions evolution at small x, especially when the 
input functions are singular^ 13 ' 14 !. 

We will first discuss the phenomenological implications of our remarks. In Fig. 3, 
we exhibit the data with corresponding fits for Q 2 = 15 and 30 (GeV/c) 2 , values which 
correspond to the best-explored small- a; domain. The lack of statistics in the present data 
excludes the possibility of obtaining a precise estimate of dlnF^/dln -k For this reason, we 
will proceed as following. We consider the scale /i in eq. (8), as the only relevant parameter 
with A fixed at a reference value of 100 MeV. Then, leaving free the normalization of the 
structure function at each Q 2 , we obtain an estimate of /i = 290±60MeV. This preliminary 
analysis is an incentive for a more complete study which will be allowed by the larger 
statistics at HERA in the near future. At the very least, it confirms that either a BFKL 
or a GLAP singularity may well describe the present data. In agreement with the GRV 
phenomenological analysis^, the value of /i is situated in the non-perturbative range of 
scales for QCD. 

Some theoretical comments are now in order. 
The BFKL singularity 

This singularity is obtained as a resummation of the leading [a s ln -j;] n perturbative 
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QCD contributions coresponding to soft gluons. The measurement of proton structure 
functions in the low- a; region accessible at HERA, may provide the first experimental 
signal of the BFKL singularity predicted long ago from QCD calculations^ . However it 
is also known that this singularity contains non-perturbative contributions connected to 
the low transverse-momentum range of the gluons in the corresponding ladder diagrams. 
This is often cited as the J«t— diffusion problem^. A somewhat related problem is the 
the infrared cut-off dependence^ 15 ] of the BFKL singularity, especially when the coupling 
constant is allowed to run. Let us see what can be learned from the comparison with the 
GLAP evolution equations. 
The GLAP singularity 

This singularity is obtained as a resummation of the leading [a s lnQ 2 ] n perturbative 
QCD contributions corresponding to the resummation of collinear gluons. More precisely, 
the use of this singularity at small x is nothing but the double logarithm approximation 
(DLA). The DLA approximation is not expected to be dominant in the present range of 
HERA experiments Thus it came as a surprise that the GRV parametrization^ based on a 
fit of previous experimental results using the GLAP evolution equations over a large range 
of scales yields a good prediction of HERA data^ 1 ' 2 !. We now see by our method using 
dlmF/dln-j: that it is not a mere coincidence but a consequence of the striking similarity 
between GLAP and BFKL singular behaviour in the experimental small- a; range. 

We obtained our results for evolution equations including a non-perturbative range, 
fi> < Q < Qo- It is worthwhile to compare them to those obtained some time ago^ 13 ! in the 
deep-perturbative regime, Q ^> Qo. Starting with the same input at Qo, the BFKL and 
GLAP evolution of the structure functions is found to be quite similar. The theoretical 
origin of this behaviour is the approximate equality of gluon anomalous dimensions when 
the expansion parameter ^ is small^ 13 !. In our analysis, the expansion parameter near 

j ' = 1 is j5±, see eq. (1), remaining of order 1 in the domain of phenomenological interest. 
Our conclusion is that the end-point of the singularity a is compatible with a genuine 
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BFKL singularity, while the type of singularity is modified. This can be seen from Fig. (2), 
where the parameter b evolves from a value not far from the predicted one (b ~ — |), 
towards positive values at larger ^. 

One may object to the validity of the GLAP formalism in a region which includes 
non-perturbative contributions^ 16 ^. But, as we have already mentioned, non-perturbative 
contributions also have to be taken into account in the BFKL formalism for deep inelastic 
structure functions^ . There may well be less difference between the two approaches at 
small x than expected. Let us discuss this point a little further. If this similarity is more 
than a coincidence, one obtains a new information on the non-perturbative region of the 
kT~ diffusion of the BFKL gluon radiation. This suggests that the dominant ln-j:,lnA;x 
ranges for gluons in the ladder diagrams^ contributing to F s (assuming the extension 
of the BFKL evolution equation to the non-perturbative domain) could be in one-to-one 
correspondence with those in the GLAP equations with an evolution from a low /i 2 scale. 
The hadronic final states in Deep-Inelastic Scattering could be a sensitive test of this 
conjecture. 

As a final remark concerning the Q 2 dependence, there seems to be a slight theoretical 
contradiction between the predictions of the BFKL and GLAP evolution equations at small 
x. Indeed, the latter can be proven to be valid if the input structure functions are singular 
enough^ 13 ' 14 ] while the former predicts a different (^-evolution pattern^ 12 !. We have found 
a specific Q 2 evolution of the BFKL parameters a(C),b(^) which remains to be studied 
theoretically and phenomenologically in order to clarify this point. 
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FIGURE CAPTIONS 



Figure 1: Comparison of BFKL and GLAP singularities for dlmF/dln-j: 

The predictions for dlmF/dln-j: in the small- a; range are displayed for the GLAP singu- 
larity (dotted line) and the BFKL singularity (continuous line). The small- a; range is lim- 
ited from above (excluding the dashed lines) by the range of validity of the j— plane singu- 
larity dominance (see text). The curves are displayed for three typical values £ = (.1, .2, .3) 
of the GLAP evolution parameter. The corresponding values of Q span a very large range 
between 1.7 and 810 5 GeV(!) with the chosen parameters (see formula (8) in the text). 

Figure 2: Correspondence between BFKL and GLAP singularity parameters 

The a and b parameters of the BFKL singularity, eq.(l), are displayed as functions of 
£, the GLAP evolution parameter (full black stripes). The width corresponds to systematic 
variations on the fitting range, see text. The grey region corresponds to a simple power-law 
fit (b = 0). 

Figure 3: Phenomenological comparison 

The data for the singlet structure function at small x from H\ (Q 2 = 15, 30 GeV 2 
from ref. [1]) and Zeus ( Q 2 = 30 GeV 2 from ref. [2]) are displayed in a InFs/ln 1 plot. 
The continuous line corresponds to a fit with /i = 290MeV, see formula (8). (The absolute 
normalization is fixed independently for each set of data, since we want to parametrize only 
the derivative dlmF/dln-j:). An estimate of the admissible values give \i = 290 ± 60MeV, 
which correspond to £ = .135 ± .015 ( Q 2 = 15 GeV 2 ) and to £ = .15 ± .02 ( Q 2 = 30 
GeV 2 ), all other parameters being fixed as in the text. 
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